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Drude weight of optical conductivity is calculated at zero temperature by exact diagonalization 
for the two-dimensional t-J model with the two-particle term, W . For the ordinary t-J model 
with W=Q, the scaling of the Drude weight D oc 5^ for small doping concentration 5 is obtained, 
which indicates anomalous dynamic exponent z—i of the Mott transition. When W is switched 
on, the dynamic exponent recovers its conventional value z=2. This corresponds to an incoherent- 
to-coherent transition associated with the switching of the two-particle transfer. 
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High-Tc superconductors (HTSC) are a typical exam- 
ple of strongly correlated electron systems, and show 
many unusual physical properties which challenge, the 
picture of the standard Fermi liquid theory.ErBQ'cP In 
particular, transport properties are among the most pe- 
culiar ones. Electric resistivity exhibits a linear temper- 
ature dependence p{T) ~ T near the optimal doping, 
and the frequency dependence of optical conductivity 
is long-tailed, cr(w) ~ uj~^. These are not compatible 
with the standard Fermi liquid theory, which predicts 
p{T) ^ and (j{lo) ^ and indicate that electric 

transport in HTSC is anomalously incoherent. Similar 
incoherent transport properties have also been observed 
in many other strongly correlated systems including a 
number of transition-metal oxidestP and organic conduc- 
torSjtP and therefore it may be natural to ascribe the 
incoherent character to strong correlation effects. 

Concerning the mechanism of high-Tc superconductiv- 
ity, many experiments indicate that the strong electron- 
electron correlation may play an essential role, and vari- 
ous theoretical approaches have been proposed with, 
phasizing different aspects of the correlation effects.llr 
As discussed above, incoherent motion of electrons is a 
notable ingredient of the strong correlation effects, and 
one possible scenario is that the incoherence in the nor- 
mal phase may drive the superconductivity so that the 
system can acquire an additional energy due to coher- 
ent matipn which is recovered in the superconducting 
phase.iB 

This idea has recently been explored in .terms of the 
quantum criticality of the Mott transitionjj'' where sin- 
gularity of two quantities has been numerically investi- 
gated- pOne is the doping dependence of compressibil- 
ity kJU'I^ and the other is the localization length of the 
single-electron Green's function, ^j, pyith varying chemi- 
cal potential in the insulator phaseu-' Both results sup- 
port the scaling hypothesis of the Mott transition with 
dynamic exponent z=4, larger than the conventional 
value z=2. If this is the case, carrier transport is ex- 
pected to be largely incoherent. However, this has not 
been examined directly. It is imperative, for the direct 



experimental relevance of the above idea, to clarify the 
incoherence of transport properties. Moreover the singu- 
lar behavior of compressibility could allow an alternative 
interpretation, the spin density wave (SDW) or equiva- 
lently the nested Fermi surface. 

In this paper, we investigate coherence of electronic 
transport in strongly correlated electron systems by nu- 
merical calculations for the t-J model. First, we cal- 
culate the doping dependence of the Drude weight at 
zero temperature. The result shows the failure of the 
SDW picture and supports the scaling theory. We also 
determine the dynamic exponent z to quantify how co- 
herence is suppressed on approaching the Mott transi- 
tion. It turns out that the spatial dimensionality is cru- 
cial and the two-dimensional (2D) t-J model is much 
more incoherent than the ID case. Secondly, we ex- 
plore whether the incoherent nature in the 2D system 
could be changed to evolve coherence by including ad- 
ditionaJ-elpitron-electron interactions. Recently, Assaad 
et aZ.t3llII'0 found an insulator-superconductor transi- 
tion in the 2D half-filled t-U-W model, i.e., the Hubbard 
model with additional two-electron processes. Supple- 
menting the t-J model with a similar term, we calculate 
the Drude weight and determine z to see if coherence 
changes. |_. 

The Hamiltonian to investigate is the 2D t-J modelEj^ 
and its variant generalized by including the square of 
local kinetic energy. For later convenience, we include an 
Aharanov-Bohm (AB) flux, ^, and then the Hamiltonian 
using standard notation is given as 

Ti = Htj + T-Lw, (1) 
Uu = -tY^ P(x.,,cLc,.+H.c.)P+J^ 5, ■ S„{2) 



Uw = ~wY,P^{xj,,+sc],c,+s.+^.c.)\%. (3) 

Here V is the Gutzwillcr projection operator excluding 
double occupancy at every site, and the sum should 
be taken over nearest-neighbor sites. The phase factor 
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is given by Xi,j = exp[— i0 • (r^— rj)/a(/)o], where (fiQ is 
flux quantum divided by 27r, a is the lattice constant, 
and we shall set a=l, h=l, and the velocity of light c=l 
throughout this paper. It should be noted that <p is 
constant, corresponding to a pure gauge, and therefore 
the magnetic field is zero everywhere. _ _ _ 

One may derive the main part of the VF-terrr£3'liil'E3' 
by the strong coupling expansion of the Hubbard model, 
and then W oc P/U with U being the on-site Coulomb 
repulsion. It is also useful to rewrite the Il^-term by 
expanding the square in eq. (3) taking account of the 
projection operators. The result is 



Vi2c] 



J<7 



(4) 



where C, is the number of nearest-neighbor sites (C = 4 for 
a square lattice) and N^. is the number of electrons. We 
note that as shown iri_eq. (4), the M^-term contains so- 
called three-site termsEir and also renormalizes the value 
of J in TLtj- 

Optical conductivity at r=0 may be defined via the 
Kubo formula asO 



V LO- 



n] 



2{E^-Eo)\{a\J, 



(5) 

where — e is the electron charge, ^ denotes the spatial 
index, V is the number of sites, and r/ is an adiabatic 
constant. |q;)'s are eigenstates of the Hamiltonian with 
eigenenergy where |0) denotes the ground state, and 

We may separate the Drude part from the formula (|^) 
and the remaining is the regular part: 



V 



The Drude weight is simply given byi) 



D, 



l(«l^;|0)P _ 1 d^Eoict>) 



Ea—Eo 



(7) 



(8) 



In addition to the Drude weight, another important 
quantity is the total weight of the optical conductivity. 
It is straightforward to rcwritaihe Kubo formula (||) and 
obtain the following sum rulal3^ 



(9) 



It is noted that if W—0, (F^f^) reduces to the kinetic 
energy along the /i-direction, while it does not if W^O 
or if long-range electron hoppings are introduced. 

To calculate the Drude weight numerically, we have 
first obtained the ground state of 7i by the Lanczos 
method, and then employed the first relation in eq. (||). 
For the latter step, we have used the Lanczos method 



again. The result agrees well with the value calculated 
via the second relation in eq. (^), and the relative error 
is typically as small as 10~^. The lattice used in this 
work is a 4x4-site cluster with periodic boundary con- 
ditions in both directions, ^==(0,0) and the total wave 
number K = (0,0), unless explicitly mentioned. Due to 
the symmetry of the cluster, the conductivity is isotropic: 
CTp^ = a, Df_t = D, and {-Ff_,^,) = {-F)Sf_,^'. 

First we show the results for the 4x4-site t-J model 
(i.e., VF=0) at various J /Vs. The Drude weight D and 
the total weight {—F) are plotted as a function of elec- 
tron density Ue = Ne/V in Fig. 1(a). The total weight 
{—F) has a maximum at ne~0.5, and it is nearly sym- 
metric around it. In contrast, the Drude weight is con- 
siderably reduced from (—F) in the region of ne>0.5. 

Similar calculatimis were performed for a 2D t-J model 
by Dagotto et alx3 to study the dependence on hole 
doping ^=1— Tie, and they concluded ~ (5 as 5 — > 0. 
We have reexamined the doping dependence and found 
a different scaling. To this end, it is convenient to cal- 
culate first the (5-dependence of the ratio of D to {~F) 
and then transform the result to that of D. Fig. 1(b) 
shows D/{—F) calculated from the data in Fig. 1(a), as 
a function of electron density. It is clear that this ratio 
becomes smaller with decreasing hole doping, indicating 
D/{—F) ~ (5 as ^ ^ 0. The kinetic energy {—F) contains 
not only the coherent part of electron motion but also the 
incoherent part. Therefore it is natural that (— F) ~ S, 
and the data shown in Fig. 1(a) confirms this scaling. 
Combining these (5-dependences, we obtain D ^ 5^ . 

The result jfor the ID t-J model is very different from 
the 2D case. til' We plot in Fig. 2 the electron-density 
dependence of D and {—F) for various system sizes at 
J/<=0.3. Difference between the two weights, i.e. the 
integrated weight of a'^'^^{ijj)^ is nearly zero at this scale. 
Therefore, not only {—F) ^ (5, but also D 5 as 5 ^ in 
the ID case. Ihis result is consistent with the Luttinger 
liquid theory.ta-' It predicts D — 2vcKc, where Vc is the 
charge velocity and asymptotically Vc ^ S. Kc is the 
Tomonaga-Luttinger liquid parameter, and Kc ~ ^ in 
the strong correlation limit, leading to D ^ S. 

The different (5-dependence of D between ID and 2D 
cases results in different values of the dynamic exponent 
z. The dynamic exponent z is defined by the divergence 
of correlation time r and correlation length ^ on ap- 
proaching the critical point: j: ~ Scaling theory for 
the Mott transition by ImadaB predicts D ~ 5i+(^-2)/<i^ 
where d is the spatial dimension. As shown above, 
D ^ 5^ and 5^ for 2D and ID t-J model, respectively, 
and this implies z=4 for the 2D and z=2 for the ID 
cases. The dynamic exponent is a measure for charac- 
terizing the coherence of the system: the larger it is, the 
faster the coherent part of conductivity vanishes at the 
Mott transition, indicating that the system is more inco- 
herent. Thus, coherence is very sensitive to the spatial 
dimensionality of the system, and the 2D t-J model is 
more incoherent than the ID model. 

The present result for z is consistent with other nu- 
merical studies, if one is based on the scaling theory. 
One is the compressibility of the 2D Hubbard modeljj' 
where k 6^^ near half filling. A similar behavior was 
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also found for the 2D t-J model by KohnoB Accord- 
ing to the scaling theory, these results indicate z—A and 
agree with the present study. Another support is the 
chemical-potential dependence of the localization length 
^/.Q^ It is found that ^; |/i — /id"'' with and 
combining with the relation vz—1, this also gives z—A. 
An alternative explanation for the (5-dependence of k 
may be given by the SDW picture, which gives the en- 
ergy E{5) ~ l\5+A5^ in the 2D case, with A being the 
momentum-independent quasiparticle gap. This implies 
K = [(1—5)^ E" {5)]~^ ~ , consistent with the numer- 
ical results. In the SDW picture, however, the kinetic 
energy is carried by coherent quasiparticles, and there- 
fore both the Drude weight and total weight should have 
a (5^-term, which is not the case as shown in Fig. 1. 

Next, we study the stability of the incoherent phase 
with z=4 found for the 2D t-J model and study what 
type of extra processes could drive the evolution of the 
coherence. A similar idea has been £xn]nied by Assaad 
et al. for the 2D Hubbard model,0'BEi where Hw 
of eq. (4) without the projection operator V was intro- 
duced. 

We show the result for the 2D t-J-W model in 
Fig. 3(a). The electron-density dependence of the Drude 
weight and total weight is plotted for J/t=0.3 and vari- 
ous values of W . For iVe=14 and 12, the ground state is 
obtained for the AB flux 0=(f, f ), and this point will 
be discussed later. With increasing W, the Drude weight 
and the total weight both increase. This is because the 
W^-term introduces additional processes of transfer in- 
volving eight second- neighbor sites, which yields a gain 
of the corresponding energy described in eq. (4). 

An important effect of the ly-term is that it affects 
the ratio D/{—F) noticeably as shown in Fig. 3(b). In 
particular, the ratio increases drastically near half- filling, 
and seemingly converges to a finite value as ^ — > 0. This 
is in contrast to the case for M^=0, where the ratio van- 
ishes linearly with 5. If this is the case, it means that 
the scaling of the Drude weight is now D ^ 5^ as in the 
ID t-J model, distinct from £> - (5^ at W^Q. There- 
fore, the dynamic exponent is correspondingly modified 
to z=2, the conventional value as in uncorrelated electron 
systems, and we can infer that the system is stabilized 
by gaining an energy of coherent motion driven by the 
M^-term. A similar change in z has also been observed 
for the 2D t-U-W modeljl3) where the exponent for the 
single-electron localization length changes from to 



It is interesting that the Drude weight calculated at 
0=(O,O) is negative at A^e=14 for W/t > 0.3. This means 
that a finite AB flux is induced spontaneously. To see 
this in more detail, we have calculated the ^-dependence 
of the lowest eigenenergy for all the sixteen wave vectors 
K. For the system size of 4x4 sites, the system with 
</)=(^, ^) (p,g: integer and, is in units of (/iq) has the 
same energy spectrum as that with 0=(O, 0). Therefore 
it is sufficient to vary ^ in a trianglular region defined 
by the three points, (0, 0), (f , 0), and (f , f ). We have 
kept track of the lowest eigenenergy for all the K^s along 
the edges of the trianglular (/)-region, and the result is 
shown in Fig. 4. The eigenstates have the lowest energy 



at ^=(f,f) and K=(-f,-f), when W/t > 0.05. A 
finite AB flux is also found for iVe=12 when W/t > 0.1, 
and the value is again (j)=[j, ^), but now at K={'k, tt). 

This spontaneous AB flux may be an evidence of en- 
hancement of d2.2_j^2-wave superconducting correlation 
due to the M^-term. Actually, such a superconducting 
phase has been observed in the t-U-W model. O'O Re- 
call that the flnite flux shifts one-particle momentum as 
k—>k^<p, and no electrons then have the renormalized 
momentum k' — (±-|,±-|). Therefore all electrons on 
the tight-binding "Fermi surface" can utilize a pairing 
potential with d^-2_y2 symmetry. This also suggests the 
importance of strong momentum dependence of charge 
excitations, and that incoherence around (±7r, 0) and 
(0, ±7r) rather than (±^,±^) is responsible for the re- 
trieval of the coherence and the pairing. We will discuss 
this point in more detail elsewhere. 

In this paper, we have numerically studied coherence 
character in transport properties of strongly correlated 
electron systems. We have calculated the doping depen- 
dence of the Drude and total weights of optical conduc- 
tivity for the 2D and ID t-J model at zero tempera- 
ture. The result shows different scalings near half filling, 
D and for the 2D and ID case, respectively. The 
dynamic exponent of the Mott transition as (5 —> is 
determined as z=A for the 2D case while z—2 for the ID 
case. Incoherent transport properties in the 2D system 
may be attributed to this anomalous dynamic exponent. 
We have also examined the effects of the extra electron- 
electron interaction, Ti.w, in two dimensions. It turns 
out that the dynamic exponent changes to 2=2 and co- 
herence is retrieved. 
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Fig. 1. (a) Drude and total weights of the optical conductivity of 
the 2D t-J model with 4x4 sites at r=0. (b) The ratio of the 
two weights D/{—F). 



Fig. 2. Drude and total weights of the optical conductivity of the 
ID t-J model at T=0 for J/t=0.3. 



Fig. 3. (a) Drude and total weights of the optical conductivity 
of the 2D t-J-W model with 4x4 sites at T=0 for J/t=0.3. D: 
open symbols, and (—F): solid symbols. For ne=^ and y|, the 
ground state is calculated at (t>={^, ^). (b) The ratio of the two 
weights D/{-F). (c) The VK-dependence of the ratio D/{-F) 
for ne=|| at J/t=0.3. 



Fig. 4. The AB-flux dependence of the lowest eigenenergy for all 
wave vectors for the 2D t-J-W model, 14 electrons in 4x4 sites, 
and J/t=0.3, W/t=1.0. 
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